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It is shown that the mode-coupUng equations for the strong-couphng hmit 



S^ . of the KPZ equation have a solution for d > 4 such that the dynamic ex- 

ponent z is 2 (with possible logarithmic corrections) and that there is a 
delta function term in the height correlation function {h(k,uj)h*(k,uj)) = 
{A/ k^'^^^)6{uj / k^) where the amplitude A vanishes as d ^ 4. The delta 
function term implies that some features of the growing surface /i(x, t) will 
persist to all times, as in a glassy state. 
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A simple non-linear Langevin equation has been proposed by Kardar, Parisi and Zhang 
|jl|] and is now widely accepted as describing the macroscopic properties of a wide variety of 
growth processes, such as the Eden model, growth by ballistic deposition and the growth of 
an interface in a random medium 0. This equation is also related to other seemingly dis- 
parate problems such as randomly-stirred fluids |^ (Burgers equation), dissipative transport 
in the driven-diffusion equation 0], the directed polymer problem in a random potential § 
and the behavior of flux lines in superconductors ^ . Because of its ubiquity, any advance 
in understanding the KPZ equation is likely to have wide significance in both the fields of 
nonequilibrium dynamics and disordered systems. 

The KPZ equation for a stochastically growing interface is: 

^^ = -^'h+ ^(Vhr + r^(x,t). (1) 

It describes the large-distance, long-time dynamics of the growth process specified by a 
single-valued height /i(x, t) (i.e. one with no overhangs or voids) on a ci-dimensional sub- 
strate, X G W^. This equation reflects the competition between the surface tension smooth- 
ing forces, z/V^/i, the tendency for growth to occur preferentially in the direction of the local 
normal to the surface, represented by the nonlinear term in Eq. (P and the Langevin noise 
term rj which is added to model the stochastic nature of this growth process. The noise has 
zero mean and is Gaussian such that 

(r7(x,t)r/(x',t')) = 2D6\^-^')6{t-t'), (2) 

where D specifies the noise amplitude. 

The objective is to characterize the form of the surface. Commonly studied are the 
correlation function 

C{k,uj) = {h{k,uj)h*{k,uj)) (3) 

and the response function: 

^ '"^^ ~ 6^{k~k')6{uj-uj')\dv{k',uj')/- ^ ' 
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The correlation and response function take the scahng forms: 



c(k.-) = pjki"(Fl (^) 



''<''•'"> = FnF'- ('') 



For d > 2, there are two distinct regimes. There is a weak-couphng regime, for A < Ac, 
where the nonhnear term is irrelevant and z = 2 and x = (2 — (i)/2. For A > Ac, the 
non-linear term is relevant and the scaling relation x + 2; = 2 follows from the invariance of 
Eq. (|I]) to an infinitesimal tilting of the surface /i— s>/i + v-x, x— i>x — Avt 0. There is 
thus only one independent exponent to be determined in the strong-coupling regime, (which 
is the only one we shall consider here). 

Because there are no obvious small parameters to describe the strong-coupling regime, 
most studies of it have been numerical [|^ . The most recent numerically determined values of 
the dynamic exponent z seem to lie between the Wolf-Kertesz |Q conjecture xjz = l/{2d+l) 
and that of Kim and Kosterlitz [^ x/^ = ^/{d + 2) when d < 4. For d > 5, xj z is still 
apparently non-zero but lies below the values predicted by both conjectures. Thus there is 
some very weak numerical evidence that the upper critical dimension t/c beyond which z = 2 
and X = is 4. In our studies we find (ic = 4 and we beheve that for d > 4 the apparent 
non-zero values of x arise from logarithmic factors masquerading as small powers. 

We take a non-perturbative approach to the strong-coupling regime called mode-coupling 
theory |^ ; in it one retains in the diagrammatic expansion for C and G only diagrams which 
do not renormalize the three-point vertex A. This procedure leads to the coupled equations: 

G-^(k,c.) = Go-^(k,.;) + A2|^|^[q-(k-q)][q-k]G(k-q,a;-fi)C(q,fi) (7) 

\2 „ JO „ Jd„ 

C(k,^) = Co(k,^) + y|G(k,^)ry— y^[q-(k-q)]^C(k-q,.;-l])C(q,fi), (8) 

where G^i^.uj) = {uk"^ — iu)~^ is the bare response function and Co(k, cu) = 2D|G(k, u;)p. 
Some of us [jlO| have recently shown that the mode-coupling equations arise from the large- A^ 



limit of a generalization of the KPZ equation to an A^-component model. In principle, this 
might allow a systematic expansion in 1/A^, by which one could go systematically beyond 
mode-coupling equations towards a solution to the full problem. 

In the strong-coupling limit, the scaling functions n{x) and g{x), where x = uj/k^, satisfy 
the equations 

g^'^ix) = -ix + PiIi{x), (9) 

n{x) = ipi|(7(x)p/2(x), (10) 



where Pi = X^/{2'^r{^)'K'^+^/^). The integrals h and h are given by 

/•TT /'H-OO PCX) 

h{x)= dOmi'^-'^e dy dq cos 9 {cos 9 - q) 

Jo J-oo Jo 



f^ ( x-q'y \ 



/•TT /' + 00 roo 

l2{x)= d9sm'^-^9 dy dq{cos9 - qf 

Jo J-oo Jo 



x^"(^)"fe). (12) 



1/2 



where r = f 1 + g^ — 2g cos ^j and A = d + A — z. 

In the strong-coupling limit, the bare term in D in Eq. (^ can be dropped, as can the 
term uk"^ in the bare propagator Go- Eqs. (|9|-p!2D are valid for the limit uj —>■ 0, k^ —>■ with 
u/k^ fixed. Notice that provided 2; < 2 no cutoff is needed. 

The numerical solution of the mode-coupling equations in Eqs. (^ -|1^) presents 



formidable problems. Recently Tu |TT| has attempted such a numerical solution, but the 
dependence which he obtained for z on d (first increasing from the exact value oi z = 3/2 
in d = 1 then decreasing at larger d) is so strange that we suspect that his solution cannot 
be accurate. We suspect from our own attempts at finding a direct numerical solution that 
problems can arise from the integrable singularities in eqs. (|1T]) and (p^. However, if we 
first assume that dc = 4: so that for d > dc, z = 2 and that n{x) = A6{x) then progress is 
possible. (We shall later confirm that such assumptions are consistent.) 



li z = 2, it is now no longer possible to drop uk"^ from the bare propagator. Further- 
more, /i(x) ceases to be well-defined as the final momentum integral diverges logarithmically 
without a cutoff A. Thus even ii z = 2 for d = A, the scaling of u is not likely to be simply 
with /c^ but with k"^ modified by some (unknown to us) power of log{A/k). We have been 
unable to make any analytic progress once cutoffs are explicitly required. Instead we shall 
study the following problem in which the bare propagator is: 



Go\k,uj) = u{z*)k' 



lUl, 



(13) 



and we shall imagine that z* is arbitrarily close to 2. With 2;* < 2 no cutoff is required. 
Setting n{x) = A6{x), eq. (|) then becomes: 



g-\x) = g-\0)~zx + Piihix)-hiO)). 



(14) 



The integrals defining the difference (Ii(x) — /i(0)) are convergent without a cutoff and can 
be calculated with z = 2. Also ^"^(0) = u{z*) + /i(0) where /i(0) = X'^Ag{0)Ti{z*, d) with 



Ti{z,d) 



(3 - 2z) r 



l)5(l-f,^-l) 



(4^)d/2+ir(| + 2 - ;z)r (^) ' 

where B{x,y) is the beta function. In addition, with n{x) = A6{x), eq. (|10D becomes: 



(15) 



A = X'A'\giO)\'T2iz*,d), 



(16) 



where 



T2{z,d) 



t(^^ + 2-2z)b{z,z-1] 



(47r)d/2+i y(I + 2-z 



- + 2{z-l){z-2) 



(17) 



Notice that Ti{z, d) is divergent as z ^ 2 reflecting the need for a cutoff in that limit. The 
solution of ([T6|), setting z* = 2 (the integrals here are convergent in that limit) is: 

(rf-4)(d-2) 4r(d/2) 



A\gm' 



p,d r(i/2)r((d-i)/2)- 

Setting g{0) = 1 (as can always be achieved by adjusting i^{z*), the equation for g{x) is: 



(18) 



9 {X) 



1 — ix — B dr r 



-B 



A-l 



X 



[ii'i^ 



(19) 



where B = 4{d — A){d — 2)/d'^. This equation is readily solved for g{x) numerically. Notice 
that as (i — > 4, g~^{x) —>■ 1 — ix and A -^ 0. 

We have therefore found an exact solution for the mode-coupling equations for d > 4 
when the bare propagator is as given by equation (^) with z* < 2. The model with 



this bare propagator (|T^) is in some sense "long-ranged" compared to the model with the 
conventional bare propagator {i.e. with z* = 2) If the value of z emerging from the mode- 
coupling approach with the new bare propagator (eq. ([T3|) ) had been less than the bare 
z*, that is, if the renormalized propagator were even longer ranged, then the new model 
(with z* < 2) and the conventional model (with z* = 2) would belong to the same (strong- 
coupling) universality class. However, this was not found. The z of the calculation is equal 
to the z* of the bare propagator. One concludes that the value of z associated with the 
conventional "short-range" propagator must then be greater than or equal to z*. But as z* 
can be taken arbitrarily close to 2, we conclude that the true value of z associated with the 
true "short-range" propagator must be 2, up to logarithmic factors. 

One can check whether the solutions for n{x) and g{x) are iteratively stable as follows. 
By writing n{x) = A6{x) + Pn+i{x) in Eq. (0), etc., one sees that: 

^^dq 



Pn+i{x) = B\g{x)\'^ 
dq 



/ X 
Pn — 
q W, 



+ 



1 q 



d+l 



dq^-{d-l) 



/ X 



+ 0{pI) 



(20) 



Under iteration we found that p„+i(x) -^ \r Pn{,x) as n — > oo, with the eigenvalue \r < 1 
(which implies stability) for d < d* < 4.76.... In fact there is a relation between the 
eigenvalue \r and the behavior of Pn{x) as a; — > 0; if Pn{x) — * D /x"' as x — > 0, then direct 
substitution into Eq. (|20D shows that: 



\R 



B 
~2 



d 



d 



d/2-l-a d/2-a 



(21) 



Xr has a minimum as a function of a. Within our hmited numerical accuracy, X^ determined 
by iteration of eq. ( |20D is equal to this minimum value. When A/j > 1, i.e. when d > d* the 
simple delta function solution is unstable. We then expect that n{x) = A6{x) +p{x) where 
p{x) is proportional to {d — d*)poo{x) and Poo{x) is the limiting form for p„(x) as n ^ oo. 

Thus we have a solution of the mode-coupling equations for d > 4, which is exact in the 
scaling limit and stable. It is a "glassy" solution, in that on Fourier transforming to (k, t) 
variables, one sees that C(k, t) (in the scaling limit) is constant in time. This is rather 
like the original definition of Edwards and Anderson [J12l of spin-glass order, i.e., the spins 



Si(t) have such order if C{t) = j^I]{Si{0)Si{t)) 7^ as t ^ cxd so that the Fourier transform 
of C{t) has a 5-function in it. However, in the present case, quenched disorder is a priori 
absent, as in a 'true' glass. Hence the KPZ equation may well be another interesting model 
where quenched disorder is 'self-generated', as recently proposed and discussed in |T^. If 



this scenario is correct, our implicit assumption that the correlation and response functions 
are time translational invariant may not be valid, and the mode coupling may have to be 
recast in a two-time formulation, [p^] . 

One might wonder if the glassy behavior is attributable to the approximations made in 
the mode-coupling equations. While our solutions of n{x) are only within the context of 
mode-coupling, it is easy to see that non-mode-coupling diagrams for C, (see Fig. 1) are 
such that if the 5-function ansatz is inserted for the correlator within the diagram, then 
each of these diagrams gives only a (5-function contribution to n{x). Moreover, explicit 
evaluation of higher-order diagrams permits a generalization of eq. ([TsD which as d — > 4"*" 
takes on the form: 



where C2 and C3 are constants. Eq. ( |2^ ) implies that provided there is a nontrivial 
solution, it will always be such that A^yl|(7(0)p ~ (d — 4). Hence we expect the upper 
critical dimension dc to be 4 even beyond the mode-coupling approximation. 

An approximate solution of Eqs. (^ and (H) has also suggested that dc ~ 3.6 [|l^ [|T5| . 



Previously Bouchaud and Georges |T^ had argued that dc > 4: based on a comparison 
with directed percolation. The existence of a finite dc is supported by a 1/d expansion 
p!7| ; in addition, a prediction that dc is 4 is contained in the functional RG calculation of 
Halpin-Healy |T^ . 



For dimensions rf < 4, we do not expect to see this (5-function, but precursors of glassy 



behavior such as very long-lived peaks in /i(x, t) are known to exist for d = 2. |T9| It would 
be valuable to do numerical studies of the scaling limit of C(k, u) for d > 4 to check the 
existence of glassy behavior. 

Finally, we speculate that given the solution for d = 4 in the mode-coupling equations 
it should be possible to construct a perturbative expansion for z in e, where e = 4 — d. So 
far, however, we have not succeeded in this aim. 

Two of us (MAM and TB) would like to thank the Newton Institute, Cambridge for its 
hospitality during the performance of this work. 
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FIGURES 

Fig. 1 Diagrams for the correlation function beyond the mode-couphng approximation. 
The hnes with the circles within them are height correlation functions. If a 5-function 
form is inserted for them, the diagrams themselves give (5-function contributions to 

n(x). 
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